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ABSTRACT 

FU Orionis (FUOR) outbursts are major optical brightening episodes in low-mass 
protostars that evidently correspond to rapid mass-accretion events in the innermost 
region of a protostellar disc. The outbursts are accompanied by strong outflows, with 
the inferred mass outflow rates reaching ~ 10% of the mass inflow rates. Shu et al. 
proposed that the outflows represent accreted disc material that is driven centrifu- 
gally from the spun-up surface layers of the protostar by the stellar magnetic field. 
This model was critiqued by Calvet et al., who argued that it cannot reproduce the 
photospheric absorption-line shifts observed in the prototype object FU Ori. Calvet 
et al. proposed that the wind is launched, instead, from the surface of the disc on 
scales of a few stellar radii by a non-stellar magnetic field. In this paper we present re- 
sults from numerical simulations of disc accretion on to a slowly rotating star with an 
aligned magnetic dipole moment that gives rise to a kilogauss-strength surface field. 
We demonstrate that, for parameters appropriate to FU Ori, such a system can develop 
a strong, collimated disc outflow of the type previously identified by Romanova et al. 
in simulations of protostars with low and moderate accretion rates. At the high accre- 
tion rate that characterizes the FUOR outburst phase, the radius r m at which the disc 
is truncated by the stellar magnetic field moves much closer to the stellar surface, but 
the basic properties of the outflow, which is launched from the vicinity of r m along 
opened-up stellar magnetic field lines, remain the same. These properties are distinct 
from those of the X-celerator (or the closely related X-wind) mechanism proposed by 
Shu et al. - in particular, the outflow is driven from the start by the magnetic pressure- 
gradient force, not centrifugally and it is more strongly collimated. We show that the 
simulated outflow can in principle account for the main observed characteristics of 
FUOR winds, including the photospheric line shifts measured in FU Ori. A detailed 
radiative-transfer calculation is, however, required to confirm the latter result. 

Key words: accretion, accretion discs - MHD - stars: formation - stars: magnetic field 
- stars: winds, outflows. 



1 INTRODUCTION 

FU Orionis systems (hereafter FUORs), named after their pro- 
totype object, are low-mass (Sun-like) protostars that un- 
dergo a rapid accretion episode in the innermost region of 
the circumstellar disc (see |Hartmann & Ken yon 1996 for a 
review). The gravitational energy released in such an event 
leads to an emission outburst (with a rise time of ~ 1 — 10 yr 
and a duration of ~ 10 — 100 yr), during which the disc 
is more luminous than the central protostar by a factor of 

* E-mail: akonigl@uchicago.edu (AK); romanova@astro.cornell.edu 
(MMR); RVLl@comell.edu (RVEL) 



r-j 10 2 — 10 3 . The inferred mass accretion rate during the out- 
burst is Mi n ^ 10 -4 Mq yr -1 , much higher than typical accre- 
tion rates during the quiescent phase. Statistical arguments, 
first advanced by Herbig (1977| ), indicate that such outbursts 
occur, on average, ten or more times during the protostel- 
lar lifetime. There is evidence that the outbursts are more 
frequent during the early (the so-called Class and Class 
I) evolutionary phases and peter out as the mass accretion 
rate declines and the protostar enters the visible (Class II, or 
Classical T Tauri) phase. The picture that has emerged from 
the observations and their interpretation is that most of the 
mass that ends up in the protostar is transferred from the disc 
during such outbursts (e.g. |Calvet, H artmann & Strom 2000 ; 
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see, however, jHerbig, Petrov & Duemmler|2003| for a differ- 
ent viewpoint) . If our understanding is correct, the FU Orionis 
phenomenon represents a key element of the star formation 
process. 

FUOR outbursts are accompanied by strong winds of 
maximum lin e-of-sight speeds > 300 kms -1 (e.g. Bastian & 



Mundt||l985|), whos e inferred mass outflow rates M 



reach ~ 0.1 Mj n (e.g.|Croswell, Hartmann & Avrett 



1987 



Cal- 



|vet, Hartmann & Kenyon| |1993| hereafter CHK93). If most 
of the mass accreted through the disc is indeed processed 
through outbursts of this type then most of the mass and mo- 
mentum ejected over the protostellar lifetime - and hence 
most of the impact that protostellar outflows may have on 
their environment (e.g. in contributing to the dynamical sup- 
port of the parent cloud against gravitational collapse and to 
the regulation of the mass inflow to the centre) - will be asso- 
ciated with these eruptions. It is also likely that the repeated, 
powerful ejections have a strong influence on the proper- 
ties and appearance of the large-scale jets that emanate from 
these protostars (e.g. |Rei purth 1990). It is therefore impor- 
tant to understand the nature and origin of these outflows. 

Neither thermal nor radiative acceleration is likely to be 
important in FUOR winds, which leaves magnetic driving as 
the most promising mechanism. This conclusion is supported 
by a Zeeman-signature least-square deconvolution measure- 
ment in the prototype object FU Ori, which was interpreted as 
indicating the presence of a ~ 1 kG poloida l magnetic field 
on scales of ~ 0.05 au ( Donati et al.||2005^ f] One possible 
scenario is that the outflows represent a centrifugally driven 
wind that is launched along magnetic field lines that thread 
the disc and are sufficiently inclined (at an angle < 60° for a 
Keplerian rotation law) to the disc surface (e.g. |Blandford &| 
Payne jl982| l. These field lines could correspond to the inter- 
stellar magnetic field that threads the natal molecular cloud 
core and is dragged in by the accretion flow, although an ori- 
gin in a disc dynamo is also conceivable. In this picture, the 
outflow need not be launched from the immediate vicinity of 
the central star. Alternatively, the wind, while still compris- 
ing material removed from the accretion disc, could be driven 
along stellar magnetic field lines. In view of the fact that the 
massive accretion flow during an FUOR outburst is expected 
to compress the stellar magnetosphere to an equatorial radius 
r m not much larger than the stellar radius R*, the wind in this 
scenario necessarily originates close to the stellar surface. One 
version of the latter scenario, proposed by |Shu et al.|l|1994[l, 
corresponds to the X-celerator model presented in Shu et ah] 
(1988). In this picture, the accretion flow spins up the outer 
layers of the star to breakup, resulting in a magnetocentrifu- 
gal wind being driven along a narrow bundle of opened-up 
field lines that emerge from an 'X-point' at the stellar equator. 
In the X-wind model described in |Shu et ah] ( |1994P , the % 
point' is associated more generally with the corotation radius 
r cor (the radius where the disc angular speed equals the stel- 
lar angular speed), which in quiescent protostars is typically 
a few stellar radii away from the stellar surface. 

CHK93 and |Hartmann & Calvetj ( |1995| ) presented argu- 
ments in favour of the 'disc field' interpretation of FUOR out- 



1 There is also indirect evidence for a magnetic field from the detec- 
tion of hard X-ray emission in the FUOR object VI 735 Cyg (Skinner 
|etal.|2009p . 



flows and against a 'stellar field' scenario. They demonstrated 
that increasingly stronger photospheric lines observed in FU 
Ori become progressively more blueshifted even as their two 
absorption components (attributed to the disc rotation) move 
closer together in wavelength, and pointed out that this is 
precisely the behaviour expected in a disc-driven wind. In ad- 
dition, based on the evidence that most of the optical contin- 
uum in FUORs is emitted by an extended disc and on the fact 
that typical stellar winds accelerate on scales comparable to 
the stellar radius, they contended that a stellar field-driven 
wind launched from the stellar surface could not reproduce 
the observations. In particular, they argued that the strong, 
but only moderately blueshifted, intermediate-strength lines 
detected in FU Ori could not originate in such a wind be- 
cause, by the time such an outflow covered a significant frac- 
tion of the continuum emission region in the disc (necessary 
for producing strong absorption), it would have already at- 
tained a high velocity (and would therefore exhibit strongly 
blueshifted lines). Another potential problem that they cited 
involves the large rotational broadening that could be ex- 
pected from a source rotating at breakup (as envisioned in 
the X-celerator picture). 

Recent axisymmetric and 3D numerical simulations of 
disc accretion on to stellar magnetospheres ( Romanova et al. 
|2009| hereafter R09) have revealed features that resemble 
the X-wind configuration proposed by |Shu et aL] l |1994P but 
that are nevertheless different on several counts. Specifically, 
it was found that such systems drive conical disc winds along 
stellar field lines that are bunched up by the accretion flow. 
However, even though these winds also originate in the inner 
disc, their launching region is not confined to the immedi- 
ate vicinity of the corotation radius, as hypothesized in the 
X-wind scenario. Furthermore, the conical winds are driven 
by the pressure gradient of the azimuthal magnetic field com- 
ponent (wound up by the differential rotation between the 
disc and the star) rather than centrifugally, and they have a 
smaller opening angle and a narrower lateral extent than X- 
winds. Interestingly, even though R09 only presented results 
for model parameters appropriate to protostars with compar- 
atively low (< 10 _6 Moyr _1 ) accretion rates, the outflows 
produced in their simulations exhibited several properties 
that could potentially mitigate the aforementioned arguments 
against stellar field-driven outflow models for FUORs. In par- 
ticular, it was found that the acceleration of a conical wind is 
more extended than in a typical (hydrodynamic) stellar out- 
flow and that its rotation speed generally decreases along the 
flow, in contrast with the initial behaviour of a centrifugally 
driven wind. These findings provide a strong motivation for 
reevaluating the viability of the 'stellar field' class of wind 
models for FUORs. 

More recent investigations, employing larger simulation 
regions and higher accretion rates, have begun to extend the 
results of R09. One notable finding of this new work, anal- 
ysed in Lii, Romanova & Lovelace (2011), is that the collima- 
tion of conical winds increases with distance from the origin 
and that they can eventually become fully collimated. In this 
paper we focus on simulations that we performed for parame- 
ters that are relevant to FUORs. Our goal is to verify that con- 
ical winds are still produced under these circumstances and 
to examine whether they could potentially account for the in- 
ferred properties of FUOR outflows. In Section [2] we provide 
analytic estimates that are used to guide our simulations and 
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we summarize our numerical scheme. In Section[3]we present 
representative results and derive the physical properties of the 
simulated flows. We discuss the implications of this study for 
FUORs in Section [4] and give our conclusions in Section [5] 
where we also outline steps toward further progress. 



2 THE MODEL 

The FUOR phenomenon has been convincingly argued to rep- 
resent an enhanced accretion episode in a protostellar accre- 
tion disc, most likely associated with an instability that arises 
from a mismatch between the mass accretion rates in the in- 
ner and outer disc regions (e.g. |Zhu, H artmann & Gammle] 
2009a| |Zhu et al.|2009b| and references therein). Accordingly, 
we set up a numerical model that simulates a non-steady disc 
accretion 'burst' on to a magnetized star. We first discuss some 
basic scaling relations that allow us to choose the appropriate 
model parameters, and then briefly describe our numerical 
model. 



2.1 Physical setup 

Our model is based on the assumption that the stellar mag- 
netic field can effectively diffuse into the inner region of the 
disc, allowing the bulk of the inflowing disc material to be 
channelled on to the stellar surface along closed magnetic 
field lines (e.g. Ghosh & Lamb 1979a b) and the remainder 
to be expelled in an outflow along opened-up field lines (e.g. 
Lovelace, Romanova & Bisnovatyi-Kogan|| 1 995 1 [Goodson & 
Wingl ee|l999 1. The disc truncation (or 'magnetospheric') ra- 
dius r m corresponds to the location in the disc where the 
torque exerted on the disc plasma by the stellar magnetic field 
becomes large enough to brake the disc Keplerian rotation 
and enforce corotation with the star. For an aligned dipolar 
field, it is given by 



fci(GM» 



-1/7 



,,-2/7 4/7 
M in Ml l 



(1) 



where M* is the stellar mass and is its magnetic dipole 
moment, G is the gravitational constant and the mass accre- 
tion rate is measured at r m ( |Ghosh &~ Lamb 1979a). Under 
stationary conditions, the numerical factor ki is estimated to 
be ~ 0.5 l |Ghosh & Lamb|1979b{[Lon g, Roma nova & Lovelace] 
2005 ). When r m is close to R t , as in the FUOR case, higher 
order magnetic moments fi n , which produce magnetic field 
amplitudes B n ~ /i n /r n+2 , can also be expected to play a 
role. In this case, equation ([l| generalizes to 



= k n fi 
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-2/(4n+3) 



(GM.) 



-l/(4n+3) 



(2) 



(Long, Romanova & Lam b|20lT| ), where, in particular, n = 1, 
2 and 3 correspond, respectively, to the dipole, quadrupole 
and octupole field components. For a given dipole compo- 
nent, the incorporation of additional multipole components 
will tend to increase the value of the disc truncation radius 
over the estimate |l]l . However, for the sake of simplicity, we 
restrict the discussion in the rest of this paper to a purely 
dipolar field. 

In choosing our model parameters, we adopt as fidu- 
cial values the physical parameters inferred from observa- 
tio ns of FU Ori. In particular, based on the results given 
in Izhu et all J2007t, we take M* = 0.3 M© and M in = 



2.4 x 10~ 4 Mq yr _1 . These authors also estimate, from spec- 
tral modeling, that the inner radius r- m of the FU Ori disc is 
5Rq. We identify this radius with r m , which allows us, using 
equation ([!]), to infer the value of n n /i?*: 

■m( M. V 1/12 ( B, ^ 1/3 
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0.3 Mq 
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Mi, 



2.4 x 10- 4 Mq yr- 



-1/6 
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where the normalization of the equatorial surface magnetic 
field J5» — fii / R'i is consistent with typical values inferred in 
quiescent Class-I Oohns-Krull et al.||2009| ) and Class-II (e.g. 
|Johns-Krull||2007|) protostars as well as with the results re- 



ported by Donati et al.| ([2005 ) for the poloidal field near FU 
Ori j^] The contribution of higher order multipole field compo- 
nents, which could become important near the stellar surface, 
would have the effect of increasing this ratio [^] 

Previous accretion disc models of FUORs have generally 
ignored the role of the magnetic field in truncating the disc 
and therefore identified n n , the inner radius of the disc, with 
the stellar radius (e.g. |Zhu et al.|2008| l. The value R* = 5 Rq 
inferred in this way is measurably higher than typical values 
for low-mass protostars (< 2.5 Rq), and several explanations 
have been advanced to account for the difference. (A brief 
summary of this issue is given by |Zhu et al.|2007| who favour 
an interpretation that attributes the larger radius to stellar 
expansion brought about by the deposition of heat produced 
by the accreting gas.) Since r^fR* is generally > 1 in the 
magnetic accretion model (for example, it is 1.4 in the rep- 
resentative simulation presented in this paper), the inferred 
value of R t is lower in this case, which reduces the implied 
difference from the radii of quiescent protostars. 



2.2 Numerical setup 

We employ the same numerical model as the one described 
in section 2 of R09, and the reader is referred to that paper 
for further details. As explained in section 2.3.1 of R09, the 
high-density gas that comprises the disc material enters the 
simulation region through the disc boundary only after the 
computation commences, and it subsequently flows inward 
on account of its viscosity. This numerical setup is thus natu- 
rally suited for modelling the evolution of an accretion 'burst', 
which is the focus of the present work. 

Although we performed simulations for a variety of 
model parameters, we present only one representative case 
in this paper. We use the same parameters as in the reference 
simulation shown in R09, except that we reduce the reference 
radius Ro from Ro = 2R, to Ro = J?» and change the outer 
radius of the computational domain from i? ou t = 16 Ro to 



2 Note that the actual magnetic field strength at n n is larger than 
that of the corresponding unperturbed dipole field on account of the 
compressional amplification of the stellar field by the accretion flow. 

3 Adopting the same fiducial parameter values as in equation {3} 
and setting, for definiteness, &2 = &3 = 1-0, we deduce, using equa- 
tion {2}, that r m ,2/fi* = 1-15 and r m ,z/R* = 1.12. Note in this 
connection that recent spectropolarimetric observations of the T Tauri 
stars V2129 Oph an d B P Tau (Donati et al.|2007l|2008[ see also[Rb~l 
|manova et al.|20lT] and |Long et al.|2011| l inferred an octupolar sur- 
face field component that is stronger than the dipolar component. 
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Rout = 28i?oQThe latter change enables us to increase the 
mass accretion rate on to the central star in our simulation 
to the level inferred in FUORs. The larger size of the compu- 
tational domain also allows us to get a better handle on the 
collimation properties of conical windsj^The mass accretion 
rate on to the central object is determined from the simula- 
tions using the expression 



M simu i = M ta Mo = M in 



/'i 



p. J (GM«)V 2 Ro r/2 



(4) 



where Mo is the reference mass accretion rate and where jl 

and M are, respectively, the dimensionless magnetic moment 
and mass accretion rate parameters. We use jl — 10 as in the 

reference simulation of R09 and obtain the value of M i n from 
the numerical calculation. As described in Section [3j the fi- 
nal (quasi-steady) mass accretion rate on to the central object 

in our representative simulation corresponds to Mj n w 70. 
We also find that r^/Ro « 1.4 at that stage, which implies 
R* « 3.6 Rq (using our fiducial value for r m ). We can then 
use equation Q to infer the value of B* for our representa- 



tive model: 
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3.6 Rr. 



-5/4 
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The reference an d fid ucial parameters for our model are sum- 
marized in Table lll 6 ! 

The time evolution of the simulated system depends on 
the magnitudes of the viscosity and the magnetic diffusiv- 
ity, which are parametrized by a v and ay, respectively. R09 
(see their appendix D) found that conical winds are estab- 
lished only when a v > ad. This is consistent with the fact 
that the dragging of the magnetic field by the accretion flow, 
which in FUORs causes the field compression near the inner 
boundary of the disc, requires the magnetic Prandtl number 



Pr m = a v /ad to be > 1 (e.g. |Lubow, Papaloizou & Pringle 
|1994p . Our simulations employ the values adopted in the ref- 
erence case of R09, namely a v = 0.3 and ad = 0.1. Recent 
work on non-steady protostellar accretion-disc models (e.g. 
|Zhu et al.|2009a|b|[2010[ l has indicated that a v must be large 
enough (> 0.1) for outbursts that resemble those of FUORs to 
be produced. Our adopted value of the viscosity parameter is 
consistent with this requirement. It is noteworthy in this con- 
nection that R09 found that the formation of a robust conical 



4 For comparison, the outer radius of the outburst region in FU Ori 
was estimated by |Zhu et aL| ([2008} on the basis of infrared obser- 
vations to be r ou t ~ 0.58 au, which corresponds to ~ 35 ii* for the 
inferred value of in our reference simulation. 

5 A conical wind simulation over an even larger computational do- 
main, corresponding to Rout = 42 Ho, is presented in |Lii et aL| 
J20TTI1 



By substituting the model parameters into equation we infer 
fcx ~ 1.3. Although fci is expected to be < 1 under stationary condi- 
tions (e.g. Ghosh & Lamb 1979a}, we consider the derived value to be 
physically consistent, especially in view of the time-dependent nature 
of the simulation and the expected presence of higher order multi- 
pole field components. One could in principle obtain a lower value of 
fci by increasing the adopted value of the parameter p,. 



outflow also requires a v (and a d ) to be comparatively large 
(> 0.03). 



3 SIMULATION RESULTS 

The large-scale poloidal structure of the simulated flow is 
shown in Fig. ^ at the time when the wind has become fully 
developed. The most striking feature of the figure is its quali- 
tative similarity to fig. 3 in R09, which corresponds to a qui- 
escent protostar that accretes at a much lower (by a factor 
~ 3 x 10 -3 ) rate. This conclusion is reinforced by an in- 
spection of Fig. [2] which shows a close-up view of the re- 
gion near the star. In both cases a high-density, conical disc 
wind is launched from the vicinity of the disc truncation ra- 
dius ri n , and a lower-density, higher- velocity jet component is 
established in the interior of the cone. The main difference be- 
tween the two simulations is in the value of ri n /Ro: it is ~ 1.4 
in our simulation, as compared with ~ 2 in the R09 reference 
calculation. This difference is consistent with the expectation 
from equation pi), which indicates that this ratio scales only 
as a weak power of the mass accretion rate (MT 1 ' 6 ). Thus, 
even though the accretion flow is much more powerful in this 
case, the steep radial scaling of the magnetic pressure exerted 
by the dipolar field component (oc r~ 6 ) insures that the disc 
is still truncated at a finite radius and does not actually 'crush' 
the stellar magnetosphere (as envisioned, for example, in the 
X-celerator scenario for FUORs; see |Shu et al.|1994) . 

Although the disc in the current simulation is truncated 
very close to the stellar surface, the magnetic field structure 
in the vicinity of its inner radius is qualitatively very similar 
to the case where (r m — R t )/R t is > 1. In particular, the 
magnetic field lines that guide the conical wind and the axial 
jet are open. The opening-up of parts of the initially dipo- 
lar stellar field is a consequence of the differential rotation 
between the star, where the magnetic field is anchored, and 
the disc, into which the field lines diffuse, and is a generic 
property of magnetically linked star-disc systems (e.g 



BaUegooijen 1994; Lynden-Bell & Boily 1994| |Lovelace et al. 
1995) |Uzdensky, Konigl & Litwin|2002|. As in the lower- M in 



case considered in R09, the conical wind is driven primar- 
ily by the pressure gradient of the azimuthal magnetic field 
component generated by the differential rotation rather than 
centrifugally. Note in this connection that the profile of the 
azimuthal velocity component (see panel C in fig. 6 of R09 
as well as the right panel of Fig. [4] below) initially decreases 
along the poloidal field lines, which contrasts with the be- 
haviour of centrifugally driven winds, in which the azimuthal 
speed at first increases along a field line. The acceleration 
is quite efficient, and the conical wind reaches outflow speeds 
v p » 80 — 90kms _1 (corresponding to ~ 75 — 84% of the 
Keplerian speed at r in ) at the outer edge of the simulation re- 
gion. The highest poloidal speed observed in this simulation 
is attained further up and is associated with the lower-density 
axial flow. As seen in Fig.[T[ its value is ~ 300 kms^ 1 , which 
is consistent with the maximum line-of-sight speed measured 
in FU Ori (Bastian & Mundt 1985). However, this value may 
not be accurate since the low-density region in the vicinity of 
the axis is susceptible to numerical artifacts^ Near the base 



7 This contrasts with the robust nature of the high-velocity jets ob- 
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Table 1. Reference values (subscript '0') and fiducial values used in the representative model. See Section |Z2| and section 2.2 in R09 for further 
details. 



FU Ori 


M, (M Q ) 


0.3 


R* (R Q ) 


3.6 


Ro 


R, 


Do (cm s _1 ) 


1.3 x 10 7 


P* (days) 


7.4 


P (days) 


1.4 


P* (G) 


2.1 x 10 3 


Bo (G) 


210.0 


po (g cm -3 ) 


2.8 x 10~ lf) 


M (M yr" 1 ) 


3.4 x 10" 6 


TVo (dyne cm) 


7.0 x 10 38 




Figure 1. Poloidal matter flux pv p (in color, with the scale given at the bottom), sample poloidal magnetic field lines (yellow) and poloidal 
velocity vectors (red) for the representative conical-wind model at time T ~ 1290 d after the onset of the simulation, when the flow configuration 
in the innermost region is already fully developed. 



of the flow the wind velocity is dominated by the azimuthal 
component, which arises from the rotational motion of the 
disc in the wind-launching region and has a maximum value 
of ~ 106 km s -1 , attained at r ^ n n . (At smaller radii the 
wind azimuthal velocity decreases on account of the interac- 
tion with the magnetosphere, which rotates with the compar- 



tained in simulations of rapidly rotating stars in the 'propeller' regime 
(see R09). In the present simulation, the star is assumed to rotate rel- 
atively slowly (r cor = 3P»), which results in a comparatively weak 
axial outflow. 



atively low angular velocity of the star.) We note that, even 
though the wind is launched very close to the stellar surface 
and has a high initial rotation velocity, it does not originate 
in the stellar surface and does not require the star to rotate 
at break-up speed - which distinguishes it from the outflow 
envisioned in the X-celerator scenario ( [Shu et alTfl 988). 

As discussed in R09, the magnetic force also has a com- 
ponent directed toward the symmetry axis, which acts to colli- 
mate the wind. By using the poloidal matter flux distribution, 
R09 determined that the conical outflow in their reference 
simulation attained an opening half-angle of ~ 30° — 40°. 
From the corresponding distribution presented in Fig. [T| we 
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find that the collimation is even more efficient in the case 
that we simulate, with the outflow half-angle decreasing to 
< 10° within a radial distance (projected on the equato- 
rial plane) of < 4i?* from the stellar surface. In general, 
a magnetically driven outflow is collimated by a combina- 
tion of two effects (e.g. |Blandford &"Pa yne 1982): the mag- 
netic tension force that acts to balance the magnetic pressure- 
gradient force in the force-free sub-Alfvenic regime, and the 
hoop stress exerted by the azimuthal magnetic field compo- 
nent in the super- Alfvenic flow region]^] The difference in the 
collimation properties of the conical wind in our simulation 
and in the reference simulation of R09 can be attributed to 
the fact that a higher mass accretion rate in the disc results 
in a stronger compression of the stellar magnetic field and 
hence in a larger collimating magnetic tension force in the 
sub-Alfvenic region of the wind. The hoop-stress effect is also 
stronger in the higher- M m case on account of the compres- 
sional amplification of the field and because the differential 
rotation that twists the field lines gets stronger (for a given 
value of r cor ) as r- m is decreased. A detailed analysis of the 
collimation properties of a conical wind from a high-Af in disc 
is given in|Lii et al.| l |2011| ). 

Fig. [3] shows the evolution of the matter fluxes that are 
deposited by the accretion flow on to the stellar surface and 
in the outflow (with the mass outflow rate evaluated over 
a spherical surface far enough from the centre). It is seen 
that the mean mass accretion rate increases steadily until it 
attains a fully developed state (with A/i n corresponding to 
the value inferred in FU Ori) at a time T w 1130 d from the 
start of the simulation. This time is longer than the ~ lyr 
observed rise time of the FU Ori outburst (e.g. |Hartmann &| 
|Kenyon|[l996p , but the discrepancy is probably in large part 
just a consequence of the particular choice of initial condi- 
tions for our simulation (see Section [272] ) . In the fully devel- 
oped state, the average accretion and outflow rates are re- 
lated by Af ou t /A/i n ~ 0.13. This result is consistent with the 
observational findings in FU Ori (e.g. Cros weTTet al.|19 87). 

The rapid accretion during the FUOR outburst can be ex- 
pected to spin up the star, and it is therefore necessary to 
check whether our assumption of slow stellar rotation is self- 
consistent. We calculated the torque on the star at the end 
of our simulation from the expression N — (JVf + N m )No, 
where the reference torque No = M v R is listed in Ta- 
ble [l] and^ where die dimensionless field and matter contri- 
butions Nf and N m are, in our case, ~ 28 and ~ 23, re- 
spectively]^] By dividing the total torque calculated in this 
way, N w 3.6 x 10 40 dyne cm, into the stellar angular mo- 
mentum J* = k 2 M t R 2 A w 7.3 x 10 49 (fc 2 /0.2) gem's" 1 , 
where we assume uniform rotation with angular velocity 
fl t = (GMt/r^,) 1 ' 2 and scale the normalized radius of gy- 



8 Collimation by the ambient mass and pressure distribution, by a 
surrounding disc outflow or by an axial magnetic field anchored in 
the disc is also possible in certain cases; see discussion in R09. 

9 It is interesting to note that, whereas the field and matter contri- 
butions are comparable in this case, the field contribution strongly 
dominates in quiescent systems (see R09) . It is also worth noting that 
the angular momentum discharge to the wind is smaller than the 
torque on the star: in our fiducial case, JVf —5 and JV m « —21 at 
R = 6 R* at the end of the simulation. This is consistent with the fact 
that, as in the quiescent case (see fig. 14 in R09), most of the angular 
momentum transport occurs through the viscous stress in the disc. 



ration k 2 by its value for a polytrope of index 1.5, we infer 
a characteristic spin-up time ~ 65 (fe 2 /0.2) yr, which is of the 
order of the typical FUOR outburst time. This implies that our 
assumption of a slow rotator is only marginally consistent. 
We note, however, that a conical wind-like component may 
be present in the outflow even if this assumption is violated 
(see Section[5]). 



4 IMPLICATIONS FOR FUOR OUTFLOWS 

The simulation results presented in Section [5] indicate that 
the observed properties of FUOR outflows could in principle 
be explained in terms of the conical wind and axial jet that 
are driven from the vicinity of the stellar surface in these sys- 
tems along stellar magnetic field lines that are compressed, 
twisted, and opened up by the interaction between the ini- 
tially dipolar field component and the strong accretion flow. 
In particular, for typical values of the stellar mass, radius and 
surface magnetic field strength, and of the mass inflow rate 
at the inner edge of the circumstellar disc, our representa- 
tive simulation demonstrates that this interaction can pro- 
duce outflows whose properties (mass outflow rate, strong 
rotational velocity component near the base and possibly also 
the maximum outflow speed) are consistent with the obser- 
vations. 

As was mentioned in Section]!] CHK93 and Hart mann &| 
Calvet] ( |1995[ ) argued against a stellar magnetic field-driven 
wind being able to account for the spectral properties of the 
outflow in FU Ori. They envisioned the outflow as originating 
in the stellar surface and accelerating rapidly along strongly 
divergent field lines even as its rotation speed (which ini- 
tially has the stellar breakup value) continues to increase. In 
our picture, the absorption features modelled in the above- 
mentioned papers would arise in the conical-wind compo- 
nent, which exhibits a spatially more extended acceleration 
(along fast collimating magnetic field lines) and a lower initial 
rotation speed (that at first actually decreases along the flow) 
than the 'stellar field' outflow assumed in those papers. Al- 
though the spatial and kinematic properties of our simulated 
conical wind are distinct from the semi-analytic disc outflow 
model presented in CHK93 (which combined a hydrostatic 
disc atmosphere with a simple representation of a centrifugal 
wind), it is probably qualitatively closer to that model (which 
CHK93 and Hartmann & Calvet 1995 argued was consistent 
with the spectral data for FU Ori) than to their hypothesized 
stellar wind model. 

To demonstrate that the conical wind model could ac- 
count for the behaviour of the absorption lines measured in 
FU Ori would require a determination of the thermal struc- 
ture of the simulated flow and a calculation of the synthe- 
sized spectra of the relevant photospheric lines. An analogous 
radiative-transfer calculation, addressing the rotationally in- 
duced line variability from an accreting T Tauri star with a 
misaligned magnetic dipole, was carried out by Kurosawa, 
I Romanova^ & Harries] ( |2008| ) . While a detailed computation 
of this type is outside the scope of the present paper, we can 
obtain some indication of the potential promise of this model 
by calculating the density and velocity profiles as functions of 
distance from the mid-plane at the location of the optical con- 
tinuum emission region and comparing the results with those 
obtained in the disc wind model of CHK93. In view of the 
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Figure 2. Innermost region of the simulation domain presented in Fig. [T] showing details of the disc-star interaction and the formation of a 
conical wind. The inner radius of the disc is denoted by n n . The arrows depict poloidal velocity vectors, whereas the black lines represent sample 
poloidal magnetic field lines. The white line labelled /3 = 1 marks the surface where the gas pressure p is equal to the magnetic pressure B 2 /8n. 
The heavy dash-dotted line that starts at r = 2r m marks the trajectory along which the conical wind parameters are plotted in Fig. [4] 
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Figure 3. Time evolution of the matter discharges from the disc on to the star (top curve) and into the wind (bottom curve). The mass outflow 
rate is calculated by integrating the matter flux over a spherical surface of radius R = 6 R t . 
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clear differences between our numerical model and CHK93's 
semi-analytic model (which include the fact that the latter 
model, in contrast with our simulations, incorporates energy 
loss by radiative diffusion in the disc atmosphere), we cannot 
expect to find a full quantitative correspondence between the 
two calculations ™ However, we can look for common trends 
in the respective profiles. We take the optical emission radius 
r opt to correspond to an effective disc temperature of 5300 K, 
and we use the disc model of Zh u et aL] 1(2007), in which 
Tes(r) oc (ri n /r) 3 [l — (rh/r) 1 ' 2 ] and the maximum effective 
temperature (attained at r max « 1.36 n n ) is T max = 6420 K, 
to deduce r opt = 2.37 r in = 11.85 R Q . As noted by CHK93 , 
the disc model fits of |Kenyon, Hartmann & Hewett| ( |1988| ) 
similarly imply that roughly 60% of the optical spectrum in FU 
Ori arises from disc annuli between 1.5 n n and 3n n . CHK93 
suggested that this region could be adequately represented 
by their calculated 'disc atmosphere plus wind' structure at 
r = 2 r in , which they presented as a function of the height 
z above the mid-plane in their table 2. In view of the differ- 
ent model setups and adopted source parameters, there are 
several possible choices for the value of the wind-launching 
radius where we could compare our model results with the 
ones given in that table. For definiteness, we opt to also use 
r — 2 r in . Note, however, that this value corresponds to differ- 
ent radial distances in the two models: 10 Rq for our choice 
of parameters and 8.84 Rq for those employed by CHK93. 

The comparison between the predictions of these two 
models is presented in Fig. [4| In the CHK93 model, the 
poloidal velocity only has a vertical (z) component and the 
azimuthal velocity is approximated as being constant with 
height and set equal to the Keplerian speed at the base of 
the flow (~ 103 km s -1 for their adopted parameters). The 
curves showing v z and the mass density p in that model are 
labelled by 'CHK' and plotted as a function of z at r = 2 n n 
using the data in table 2 of CHK93. In view of the conical 
shape of our simulated wind, we plot the poloidal and az- 
imuthal velocity components as well as the density in the nu- 
merical model along a slightly inclined path (represented by 
the heavy dash-dotted straight line in Fig.[2j[^]lt is seen that, 
as expected, the values of corresponding quantities at a given 
distance from the mid-plane can be significantly different for 
the two casesj^] However, we also find that the basic veloc- 
ity and density structure of the two models is very similar. In 
both cases, the flow is rotation-dominated as it emerges from 
the disc but eventually v p comes to exceed (which, in turn, 



does not systematically increase along the flow). CHK93 esti- 
mated that the disc photosphere at r = 2 n n occurs roughly 
where the density drops to ~ 10 -9 gem -3 and v z increases 
to ~ lkms -1 . (For comparison, the sound speed in the op- 
tical emission region is ~ Skms^ 1 .) In their model, v z be- 
comes > Vtf, when the density drops to ~ 10 _11 gcm -3 . In 
our simulation, v p increases above ~ lkms -1 also roughly 
when the density drops to ~ 10 _9 gcm -3 , and v p comes to 
exceed v<j> when p decreases by another two orders of magni- 
tude. This correspondence indicates that the observed depen- 
dence of photospheric absorption line profiles in FU Ori on the 
line strength, which was successfully reproduced by the CHK3 
model (see also Hartmann & Calvet 1995 ), is consistent with 
an origin in a stellar field-driven conical wind. 

Even if more extensive and detailed calculations indicate 
that the observed behaviour of the photospheric absorption 
lines in FU Ori cannot be reproduced by a conical wind model 
(because, for example, a wind of this type that is launched 
very close to the stellar surface collimates too rapidly for its 
acceleration region to intercept a line of sight to the opti- 
cal continuum emission region for the inferred disc inclina- 
tion angle, or if the maximum predicted outflow speed is too 
low), the results of our simulations suggest that a stellar field- 
driven outflow might still be an important ingredient of a 
comprehensive model of FUORs. In particular, such an out- 
flow could still potentially account for much of the mass and 
momentum injected into the ambient medium in the course of 
an outburst and perhaps also for the highest measured veloci- 
ties in the Ha, H/3 and Na I lines (e.g. |Bastian & Mundt|1985| 
|Hartmann & Calvet||1995] ) even if another outflow compo- 
nent (in particular, a disc wind driven along non-stellar mag- 
netic field lines, which were not included in our simulations) 
gives rise to the observed photospheric lines. This is because 
a strong conical wind and a fast, low-density jet appear to 
be generic features of the disc/stellar-field interaction under 
a wide range of conditions. On the other hand, if it can be 
demonstrated that these predicted outflow components are, 
in fact, absent during FUOR outbursts, this would indicate 
that at least one of the underlying key assumptions of the 
model [e.g. that the magnetic diffusivity in the inner disc does 
not exceed the viscosity («d < a v ) but is nevertheless suffi- 
ciently large (ad > 0.03), or that the star possesses a suffi- 
ciently strong dipolar field component (p,\ > 3 x 10 37 G cm 3 )] 
is not valid, which would also enhance our physical under- 
standing of these systems. 



10 The two calculations also differ in their adopted parameters for 
FU Ori. CHK93 used M* = 0.49 Mq, r in = R, = 4.42 _R and 
M in = 1.59 x 10~ 4 A/ yr -1 , and presented results for M out = 
lO _5 M0yr -1 . Furthermore, in their spectral fits they assumed a 
disc inclination angle i = 35°, whereas more recent fits for this 
source (e.g. Zhu et al. 2007) have used i = 55°. 

11 Note that, because of the variation of the poloidal velocity with 
distance from the symmetry axis, the density distribution in the con- 
ical wind is less strongly collimated than the distribution of the 
poloidal mass flux plotted in Figs.^and|2j 

12 The vertical shift between the two moael curves in Fig. [^reflects, 
in part, the fact that our simulated disc is hotter - and therefore geo- 
metrically thicker - than the CHK93 model disc. The latter model was 
chosen to be consistent with the minimum-temperature requirement 
obtained by Clarke, Lin & Pringle ( 1990) for a viscous disc to gener- 
ate FUOR-type outbursts. The hotter disc produced in our simulation 
is clearly also consistent with this condition. 



5 CONCLUSION 

We have presented numerical simulation results that sup- 
port an interpretation of the powerful winds that accompany 
FUOR outbursts in terms of stellar magnetic field-driven disc 
outflows. In this picture, the massive accretion flow that gives 
rise to an observed burst strongly compresses the stellar mag- 
netic field lines, and the resulting magnetic stress truncates 
the accretion disc very close to the stellar surface. Some of 
the field lines diffuse into the disc and become twisted by 
the differential rotation between the disc and the star. This 
twisting, in turn, opens up the field lines, and the pressure 
gradient associated with the azimuthal magnetic field compo- 
nent along the opened field drives a moderate-velocity, dense 
conical wind that emanates from the vicinity of the trunca- 
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Figure 4. Variation of the density (left panel) and of the poloidal and azimuthal velocity components (right panel) of the conical wind along 
the ray marked by a heavy dash-dotted straight line in Fig. [2] These quantities are depicted by solid lines, starting at the location where v p first 
increases above 1 kms -1 . The dashed line continues the iu curve along the same path to lower elevations, where the disc rotates at the local 
Keplerian speed. Also shown are the variations of p and v z along a strictly vertical line segment in the CHK93 semi-analytic disc-wind model 
for FU Ori (subscript 'CHK'). In the latter model, the azimuthal velocity is assumed to be constant with height and is equal to 103 km s -1 . In 
both cases, the outflow is launched at r = 2r; n . However, because of differences in the adopted source parameters, the physical radial scale 
corresponding to this value is 10.00 Rq and 8.84 Rq, respectively, in the numerical and semi-analytic models. 



tion radius as well as a high-velocity, tenuous axial jet (whose 
properties, however, are less well determined in the simula- 
tion that we described). The magnetic field also acts to colli- 
mate these outflow components. 

The conical wind and axial jet appear to be generic fea- 
tures of the interaction between an accretion disc and a pre- 
dominantly dipolar stellar field in cases where the effective 
viscosity a v and magnetic diffusivity a d satisfy a v > a d and 
are both comparatively high (> 0.03). These features were 
originally identified in simulations of protostars with low and 
moderate accretion rates (R09). The representative simula- 
tion presented in this paper verifies that the same type of out- 
flow is produced also when the accretion rate is as high as 
~ 2.4 x 10 -4 M Q yr _1 , the value inferred in the archetypal 
object FU Ori. Our simulation implies that the disc truncation 
radius in this source, which was observationally determined 
to lie at a radius r in = 5Rq, corresponds to a distance of 
0.4 i?* w 1.4 Rq from the stellar surface, and that the sur- 
face magnetic field is ~ 2.1 kG, which is consistent with inde- 
pendent indications. The mass outflow rate in the simulated 
outburst (dominated by the conical wind) is a factor ~ 0.1 
of the mass accretion rate on to the star, and the maximum 
outflow velocity within the computational domain (attained 
in the axial jet) is ~ 300 km s -1 ; these agree well with the 
observationally inferred values for FU Ori. 

An interpretation of FUOR winds in terms of an 
accretion-disc outflow driven along stellar magnetic field lines 



was previously proposed by Shu et al. ( 1994 ) on the basis of 
the X-celerator model of |Shu et al. 119881. In this picture, 
the outflow is launched centrifugally from the surface of a 
star whose outer layers rotate at breakup speeds. This con- 
trasts with the conical-wind scenario, in which the star ro- 
tates comparatively slowly and the outflow originates at a fi- 
nite distance from the star and is driven by the B$ magnetic 
pressure gradient from the start. |Shu et aL] l |1994P general- 



ized the X-celerator model to the case where the star rotates 
below breakup, corresponding to the corotation radius r cor 
exceeding R* . However, in their generalized (X-wind) model, 
r cor still coincides with the magnetospheric radius r m , and the 
nature of the outflow from that region (the X-point) is qualita- 
tively similar to that of the X-celerator model. Our assumption 
in this paper that r m < r cor is plausible in view of the fact 
that a rapidly rotating protostar could be efficiently braked 
through a magnetic interaction with the disc during the rel- 
atively long quiescent phase (e.g. K6nig lp~991[|Ustyugova et| 
[aL]|2006] l. And while such a star would be spun up during 
the rapid accretion event comprising an FUOR outburst (see 
Section [3]), this need not result in the surface layers reaching 
breakup speeds. (Note in this connection that, even in the ab- 
sence of a large-scale magnetic field coupling the disc and the 
star, the protostellar surface layers are not expected to reach 
breakup speeds during an outburst of this type; e.g. [Popham 
|et al.|1993||T996l l. 

CHK93 and |Hartmann & Calvet| l |1995p argued that 
a stellar magnetic field-driven outflow model of the X- 
celerator type is inconsistent with the detection of moderately 
blueshifted, intermediate-strength absorption lines in FU Ori, 
which, they suggested, could be explained in terms of a disc 
outflow originating at a distance r of a few stellar radii and 
driven along magnetic field lines that are not associated with 
the star. Specifically, they showed that the observed line pro- 
files could be reproduced by a model in which gas launched 
from a Keplerian accretion disc gradually accelerates until the 
poloidal velocity component comes to exceed the azimuthal 
velocity component. In this paper we have demonstrated that 
a conical wind naturally exhibits this behaviour since the out- 
flow also starts with a predominantly azimuthal velocity com- 
ponent and eventually accelerates to v p > v^. In particu- 
lar, we showed that the density and poloidal velocity profiles 
along a ray through the conical shell that intercepts the disc at 
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the distance of the optical emission region closely match the 
corresponding profiles calculated in the disc-outflow model 
of CHK93, notwithstanding the different setups (and even 
the fiducial parameter values) employed in the two (respec- 
tively numerical and semi-analytic) models. The fact that the 
azimuthal velocity of the conical wind remains much lower 
than the breakup speed of the star and that, in contrast with 
the initial behaviour of v<j> in a centrifugally driven wind, 
it does not increase (but, rather, decreases) along the out- 
flow, circumvents another objection that CHK93 levelled at 
the X-celerator scenario. We note in this connection that an 
outflow component resembling a conical wind, as well as a 
strong axial jet component, have been found in simulations of 
the disc-magnetosphere interaction in the 'propeller' regime 
(r m > r cor ). Based on the results presented in R09, we expect 
such a flow to be more strongly influenced by the centrifugal 
force and less well collimated for given values of fj,i and M- m 
than the conical wind we considered above. However, R09 
also found in the propeller case that the magnetic force re- 
mains important in driving the wind and that the azimuthal 
speed of the wind does not increase along a field line (see 
their fig. 11). It is therefore conceivable that an outflow in 
this regime could also account for the observations, although 
this remains to be verified by an explicit simulation. 

While the results presented in this paper are highly sug- 
gestive, a more detailed calculation (involving the thermal 
and spectral properties of the outflow) is required to evalu- 
ate the contribution of a stellar field-driven disc wind to the 
absorption-line spectrum in an object like FU Ori. It will also 
be useful to carry out additional simulations in order to fur- 
ther check the dependence of the results on the adopted ini- 
tial mass and magnetic flux distributions. In particular, our 
assumption that initially there is no disc (as compared to coro- 
nal) gas in the simulation region is not realistic, and our cur- 
rent numerical setup also does not account for the possibil- 
ity that some of the stellar magnetic field may have diffused 
into the disc before the onset of the outburst (see Goodson 
& Winglee 1999). Given the comparatively high value of the 
disc inclination angle (i = 55°; |Malbet et al.|2005| l adopted 
in recent studies of FU Ori, it is conceivable that the conical 
wind model would not be able to reproduce the absorption- 
line profiles measured in this object if the source of the opti- 
cal continuum is indeed a region of size r opt > 2 r- m — 10 Rq 
in the disc. In that case a disc outflow driven along a non- 
stellar magnetic field, as proposed by CHK93, might provide 
the dominant contribution to the absorption-line spectrum of 
FU Ori, although a 'conical wind plus axial jet' outflow could 
potentially still contribute to some of the observed proper- 
ties of this object. Note, however, that if a stellar wind is also 
present (see section 5.2 of R09), it would have a decollimat- 
ing effect on the conical wind (e.g. Meliani, Cas se & Sauty| 
2006| |Fendt|2009"l ) that could increase the range of disc radii 
'covered' by the conical outflow. 

To our knowledge, photospheric line shifts such as those 
detected in FU Ori have so far not been found in any other 
FUOR. Although a lower- A/ ou t wind or some other factor 
(such as a higher projected azimuthal velocity) could have 
prevented a detection in other FUORs (see Hart mann & Cal-| 
|vet|1995] l, it would clearly be useful to be able to test compet- 
ing models also in other bursting sources. Lower-amplitude, 
repetitive photometric outbursts have been detected in EX 
Lupi and a few other T Tauri stars (e.g. |Herbig||1989| |2007| 



2008), and they have also been interpreted as enhanced mass 
accretion events. In a few of these EXOR sources there is 
evidence for an accompanying outflow, which, as discussed 
in R09, may well (at least in some cases) arise in a disc- 
magnetosphere interaction. However, these systems generally 
do not exhibit absorption-line spectra like FUORs, and there is 
also no indication that their continuum emission is dominated 
by a disc; in this regard their appearance is similar to that of 
quiescent systems. Therefore, unless other spectroscopic diag- 
nostics are identified in these sources, FUORs will remain the 
best candidates for probing the acceleration regions of proto- 
stellar disc outflows. 
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